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ABSTRACT 

Observations are revealing the ubiquity of filamentary structures in molecular clouds. As cores are 
often embedded in filaments, it is important to understand how line profiles from such systems differ 
from those of isolated cores. We perform radiative transfer calculations on a hydrodynamic simulation 
of a molecular cloud in order to model line emission from collapsing cores embedded in filaments. We 
model two optically thick lines, CS(2-1) and HCN(l-O), and one optically thin line, N2H+(l-0), from 
three embedded cores. In the hydrodynamic simulation, gas self-gravity, turbulence, and bulk flows 
create filamentary regions within which cores form. Though the filaments have large dispersions, the 
N2H+(l-0) lines indicate subsonic velocities within the cores. We find that the observed optically 
thick line profiles of CS(2-1) and HCN(l-O) vary drastically with viewing angle. In over 50% of 
viewing angles, there is no sign of a blue asymmetry, an idealised signature of infall motions in an 
isolated spherical collapsing core. Profiles which primarily trace the cores, with little contribution 
from the surrounding filament, are characterised by a systematically higher HCN(l-O) peak intensity. 
The N2H+(l-0) lines do not follow this trend. We demonstrate that red asymmetric profiles are also 
feasible in the optically thick lines, due to emission from the filament or one-sided accretion flows onto 
the core. We conclude that embedded cores may frequently undergo collapse without showing a blue 
asymmetric profile, and that observational surveys including filamentary regions may underestimate 
the number of collapsing cores if based solely on profile shapes of optically thick lines. 
Subject headings: Stars: formation, ISMistructure, kinematics and dynamics, lines and bands 



1. INTRODUCTION 

Star-forming regions are complex systems observed 
projected onto the plane of the sky. Therefore it has 
always been a difficult task to determine the true geome- 
tries and dynamics within them. One key tool for disen- 
tangling the true dynamical state of st ar-forming cores is 
the blue in f all asymmetry. Thi s effect (I Zhou et al.l 119911 : 
[Zhoul[T99l iWalker et all [1991 iMwrs et al.lll996[ ) relies 
on the fact that within a collapsing core there is both a 
centrally increasing density profile, and two points along 
a line of sight with the same velocity. In optically thick 
species (r > 1) only the foreground gas at a given veloc- 
ity is visible. This leads to an asymmetry in the double- 
peaked emission line profile as the visible blue emission 
from the far side of the core originates from higher den- 
sity gas. The blue infall asymmetry was derived under 
the assumption of spherical symmetry. 

However, observations increasingly show that this is an 
ideal rarely achieved in molecular clouds. For instance, 
I Andre et al.l (2010) show that star -forming c o res ar e fre- 
quently part of filaments, Bont emps et al.l (|2010l ) ob- 
served that massive star- formin g regions cont ain multiple 
substructures, and iTobin et^l.1 (|2010[ ) and iStutz et all 
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(|2009f ) found that asymmetries can be found on scales as 
small as 1000 AU. In this work we will revisit the asym- 
metric collapse profile using simulations of star cluster 
formation within a large molecular cloud. This paper 
examines the case of irregular cores embedded within fil- 
aments, and future work will consider the case of massive 
st ar- for n iing r egions . 

IZhoul (jl992[ ) calculate d the line p rofiles expected to 
arise from an inside-out iShd (| 19771) collapse and from 
Larson-Penston collapse (jLarsonI 119691 : iPenstonI Il969f ). 
Both models produced an asymmetric collapse profile, 
and the resulting line-widths could be used to distin- 
guish between the two collapse profiles. A follow up 
study (jZhou et al.lll993l ) showed that the line profile of 
the near-spherical Bok globule B335 matched their ex- 
pected S hu collapse profile . This result was initally con- 
firmed bv lChoi et al.l (jl995l ). but later studies have shown 
that more complex models may actually be needed (e.g. 
Wilner et al. 2000; Stutz et^l. 2008). O ther variants o f 
collapse models have been studied. Wal ker et al.l ()1994l ) 
carried out line modelling of prolate collapsing cores, not 
just sphere s, and once agai n found blue asymmetric line 
profiles. Myers et al. (1996) used a simple analytic model 
of two converging layers of gas and produced the same re- 
sult. The study of spherical core line profiles has contin- 
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ued in more recent years. In particular iRawlings fc YatesI 
(|2QQTh . iTsamis et all (|2QQ8h and Sta bler fc YenI (j2QlQD 
have shown the importance of accurately modelling the 
chemical abundances within cores, particularly since op- 
tically thick carbon bearing species, such as C O and CS, 
freez e out within the dense core centres (T afalla et al.l 

Are these models sufficiently accurate to capture the 
behaviour of realistic star-forming cores? One poten- 
tial difficulty is that for cores embedded at the end of 
a collapsing filament or other complex structure, accre- 
tion onto the core may predominantl y occu r from just 
one direction, as seen in S mith etHI (I2QTT1) . This may 
lead to either a blue or red asymmetry depending on 
which side of the core the accretion flow is located. Sim- 
ilarly, more complex morphologies and internal velocities 
within the filament and surrounding envelope could lead 
to self-absorption in the envelope obscuring the velocities 
from the core. Investigations into numerical simulations 
have demonstrated that the structures identified in syn- 
thetic maps or position-position-velocity cubes may not 
necessarily correspond to contiguous structures in three 
dimensional space, d ue to the geom e tric p ro jection along 
the line of sight fe. g. iPichardo et aIll2QQQl : lOstriker et all 
[2QQTI : IShettv et al]l2QTQ[ r 

Cores embedded within filaments are a common fea- 
ture of molecular clouds, as shown by recent Herschel ob- 
servations (Arzoumanian et al. 2011; Men'shchikov et al] 
[2Q1QI: lAn dre et al. 2010). Moreover, in Smith et alj 
(I2OIII) we specifically studied the geometry of the small- 
est scales of collapsing 'Class 0' cores from a molecu- 
lar cloud simulation, and found that even on scales of 
a ~ 1000 AU, the cores were irregular, extended struc- 
tures. The only way to understand the effects of a fila- 
mentary geometry on resulting emission line profiles is to 
directly test these scenarios with radiative transfer line 
modelling. 

Despite the many complicating factor s, blue asymmet- 
ric collapse profiles are often observed. iGregersen et al.l 
(jl997) observed a sample of 23 Class sources and 
found blue asymmetric line profiles in 9 cases, com- 
pared to 3 sources with a red asymmetry. Similarly, 
'blue' c ores have been found in s urveys by iLee et al.l 
([T999h. IGrege rsen fc Evanl (|2000h . iWu fc EvansI ^00§) 
and lFuller et al. (20051). On the other hand, while these 
surveys find potential infall candidates, the majority of 
the cores do not have a blue asymmetric profile. The 
question is, could a greater fraction of cores be collapsing 
than what is inferred by blue asymmetric line profiles? 

We focus on the line profiles arising from collapsing 
cores embedded within the irregular filaments from which 
they are formed. The layout of this paper is as follows. 
In Section [2] we describe in more detail the physical ef- 
fects which lead to a blue-peaked asymmetric line profile. 
In Section [3] we describe our methods and models, and in 
Section m we outline our results, paying particular atten- 
tion to the viewing angle of the model. In Section [5] we 
discuss how our results may aid interpretation of observa- 
tions. Finally in Section [6] we summarise our conclusions. 

2. AN ASYMMETRIC LINE PROFILE 

Figure [1] shows a schematic describing the origin of the 
blue line asymmetry in the simple case of spherical col- 
lapsing core (see also .Evans, .1999.) . The gas falls radially 
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Fig. 1. — A schematic showing the origin of a blue asymmetric 
hne profile from a spherically symmetric collapsing core. 



towards the centre, however, we only observe the veloc- 
ity component projected along the line of sight. Figure 
[1] shows how the angle (0) between the sightline and the 
core centre increases along the line of sight. Far from the 
core centre the angle is small and so the velocity com- 
ponent along the line of sight is large. However, close to 
the core centre the observed velocity component is zero. 
Conversely, the absolute magnitude of the velocity in- 
creases towards the centre of the core. For a core with 
a velocity profile shallower than Vr (x this results in 
there being two positions along the line of sight with the 
same projected velocity. 

For an optically thin line (r < 1) all the emission di- 
rected towards the observer is visible and a Gaussian pro- 
file is observed. However, for an optically thick emission 
Hne (r > 1) only the nearest point at a given velocity is 
detected. In a spherical collapsing core this means that 
on the red-shifted side of the core, the outer layers of the 
core are seen, but on the blue-shifted side, the central 
regions are seen. For an emission line with a high criti- 
cal density, there will be more emission coming from the 
core centre than from the outer layers of the core and 
consequently, the intensity is higher in the blue part of 
the line than the red part. The resulting line, therefore, 
has a central dip where the low velocity envelope ob- 
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scures the core, and two peaks from the collapsing core, 
of which the blue peak is more intense. 

The density and velocity profiles shown in the middle 
panel of Figure [T] are the key requirement to produce a 
blue asymmetric profile. The question is, are these line 
of sight profiles present in filamentary collapsing cores? 

3. METHOD 
3.1. Filamentary core models 

For our models of filamentary cores we use data from 
a giant molecular cloud (GMC) simulation that has been 
the basis of much of our recent work on star format ion. 
For details o f this s imulation see lSmith et al] (|2QQ9a[ ). In 
ISmith et all (|20TTh we found that 'Class 0' cores from 
this simulation had filamentary features on scales as 
small as a few thousand AU. In fact, in only ~ 25% of 
cases could the cores truly be considered spherical. This 
has potentially interesting consequences for core line pro- 
files, because, as discussed in the Introduction, previous 
models generally assumed spherical symmetry. 

Our original simulation used a temperature parameter- 
isation (see Smith et al. 2009b) which produced temper- 
atures of around 10 K in the dense gas (n > 10^ cm~^), 
and temperatures of up to 60 K in the most diffuse gas 
(n < 10^ cm~^). T his includes a he ating term based 
on the YSO models of iRobitaille et al.l (2006) that is ap- 
plied in the vicinity of sink particles (which represent 
the locations of star formation). In order to verify that 
the original temperature distribution of our simulation 
was not the major cause of the line profile behaviour, we 
re-ran one of our models with a constant gas tempera- 
ture of 14 K. The only noticeable change in the result 
was that the central brightness temperatures were one 
degree higher; the qualitative behaviour of the lines did 
not change. 

We select three typical regions from the GMC simu- 
lation each containing a core embedded within a larger 
filament. The cores are identified from locations of future 
sink-particle formation, and consequently our sample is 
biased as it contains only cores which truly are collaps- 
ing. We choose to select our cores in this manner so 
that we can 1) identify collapsing cores and 2) under- 
stand how reliably we can determine this when the core 
is embedded in a dense filament. The models are labeled 
Filaments A, B and C (for reference Filame nt B contains 
the fi lamentary core shown in Figure Ic of iSmith et al.l 
I2OIII ). Figure [2] shows the 850 /im dust emission calcu- 
lated from each filament. Filaments A and C are aligned 
parallel to the axis of the model grid, and each contains a 
second embedded core besides the central one. Filament 
A is a turbulent sheet. Filament B is an interface where 
fast-moving diffuse gas collides with slow-moving dense 
gas, and has a more curved geometry. Filament C is a 
head-on colliding fiow. These regions show marked simi- 
larities to recent Herschel observations f e.g. I Andre et al.l 
120101 : iMen^shchikov et aTI !2010). We discuss the dynam- 
ics and formation mechanisms of the filaments in more 
detail in Section W?2\ 

In all three cases, the embedded core contains gas with 
densities of 3 x lO^cm"^ to lO^cm"^. However when 
averaged across a Gaussian beam of half- width O.Olpc the 
maximum observed core density decreases to lO^cm"^, 
as only a small percentage of the gas in the beam has 



densities above this value. The filaments in which the 
cores are embedded have number densities of 5 x 10^ cm~^ 
to 3 X lO^cm"^. They are themselves surrounded by 
even more tenuous material, with densities of 10^ cm~^ 
or lower. 

We model the emission from a 0.4 pc box centred on 
the collapsing core. The original SPH simulation is in- 
terpolated to a grid with a cell size of 2 x 10~^ pc. For a 
core at a distance of 150 pc this corresponds to a resolu- 
tion of around 2.8 arcseconds. However, the line profiles 
are measured over a Gaussian beam with a full width 
half maximum (FWHM) of O.Olpc, which corresponds 
to 13.7 arcseconds. The models are selected at a time 
period where the central core has just started collapsing 
and a sink particle is about to form in the simulation. 
In observational parlance the cores could be considered 
as starless objects that are about to begin their 'Class 0' 
phase of core evolution. 

In order to make a comparison to observational data we 
estimate the mass of the filaments by calculating the total 
mass of gas in our grid models with a density in the range 
5xl0^cm~^ < n < 3xl0^cm~^. For Filaments A, B and 
C we obtain v alues of 6.4 Mc^ 12.3 M© and 14.7 M© 
respectively. iHacar fc Tafallal (j2Qlli ) find values in the 
range 4.8 to 11.3 for four filaments they studied 
in L1517. Our values are therefore in good agreement 
with the masses of observed filaments. 

Furthermore, we also find good agreement when com- 
paring our simulation to observed column densities in 
star-forming regions. The bottom panel of Figure [2] 
shows a column density projection from the original sim- 
ulation of each filament before radiative transfer, which 
has been smoothed as if observed by a O.Olpc width 
Gaussian beam. The column densities of our filaments 
are in the range 10^^ — 10^^ cm~^ with a typical col- 
umn density of around 4 x 10^^ cm~^, and an increase 
in colu mn density at the lo c ation of the dense embedded 
cores. IMen'shchikov etldl (|2010[ ) find column densities 
in the range 5 x 10^^ — 1.4 x 10^^ cm~^ for filaments in 
Aquila, and column densities of 3 x 10^^ — 8.6 x 10^^ cm~^ 
for filaments in Polaris. Our filaments resemble dense 
filaments in Aquila far more than the Polaris filaments. 
This is reasonable as they are extracted from a simu- 
lation undergoing rapid star formation like that seen in 
the Aquila field. The Polaris region contains no clear ex - 
amples of pre-stellar cores (Ward-Thompson et al."2010'). 
In conclusion our model filaments closest observational 
analogue is to dense star-forming filaments in clustered 
regions rather than diffuse quiescent filaments. 

3.2. Molecular line modeling 

We use the radiative transfer code RADMC-3D Q to 
carry out the line modelling of the cores. RADMC-3D 
models a variety of radiative processes including dust 
thermal emission and absorption, dust scattering, as well 
as gas atomic and molecular lines. It is this latter capa- 
bility which is of importance here. 

The chemical species which we will consider in this pa- 
per have emission lines with high critical densities and 

^ This code is publicly available with the permis- 
sion of the author, Cornelius Dullemond, at the website 
http: / / www. ita. uni-heidelberg. de / dullemond / software / ra dmc-3d/| 
A paper describing its usage is currently in preparation. 
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Fig. 2. — Top The three filaments in 850 /im dust emission. Middle The emission from Filament A for each species. Bottom A column 
density projection of the original simulation before radiative transfer, which has been smoothed as if observed by a 0.0 Ipc Gaussian beam 



as such are typically not in local thermal equilibrium 
(LTE). Consequently, a non-LTE approach must be used 
for the line transfer. We use the large velocity gradient 
(LVG) approximation ( S obolevl Il957l ) which uses local 
velocity gradients to define photon escape probabilities. 
For a full description of h ow LVG has bee n implemented 
in RADMC-3D see >Shettv et all (j2Qllal a). We also 
use a method known as 'Doppler catching' to interpo- 
late under-resolve d veloci t ies alo ng each line of sight as 
described in Shett v et all (l2Qllbl . b) 
In the cores studied here a given line of sight velocity 



may occur at more than one position and so t he assump- 
tion o f LVG is not fully satisfied. However, lOssenkopl 
(|1997l ) showed that the Sobolev approximation can be 
used in regimes where it does not strictly apply, such 
as turbulent fiows in molecular clouds, and still yield 
reasonably accurate results. In spherical homogeneous 
fiows like the idealised case shown in Figure [21 he found 
that the error had a maximum value of only 20%. More- 
over this error was mainly in the line intensity, rather 
than in the shape of the line profile. We test this fur- 
ther by calculating the LTE line profiles for our cores 
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and then comparing them to the LVG results. In LTE 
the maximum intensities are higher due to the artificially 
high level populations. However, both methods show a 
dip at the same velocity in the line profile, and so we 
can be confident that our assumption of LVG is not af- 
fecting the qualitative behaviour of the lines. We also 
tested our method using a completely spherical collaps- 
ing core for an optically thick species and confirmed that 
the expected blue asymmetric line profile was produced 
in agreement with previous studies. 

3.3. Molecular Species 

The general procedure for creating the molecular line 
profiles is the following. The filaments from our original 
simulation are mapped onto a grid and the number den- 
sity of each desired species is calculated using simplified 
prescriptions for the abundances as a function of density. 
The grids are used as inputs for RADMC-3D, which cal- 
culates the level populations of the tracer species, and 
then carries out ray-tracing to calculate the molecular 
emission. To compute the level populations of our chem- 
ical species we use line data obtained from the Leiden 
LAMDA database (Schoier et al. 2005) and assume that 
molecular hydrogen is the dominant collisional partner. 

We consider three molecular species in this study, each 
having different properties and tracing different regimes, 
as shown in Table [H The chosen line transitions all lie 
within the observable range of most ground based detec- 
tors, such as ALMA. All transitions have a high critical 
density, and hence are dense gas tracers. Consequently, 
any differences in the emitted line profiles are primarily 
due to optical depth effects. 

The isolated component of N2H+ is optically thin and 
traces the entirety of the dense gas. Molecular collisional 
rates and the the Einstein coefficients of thi s species have 
been calculated bv lDaniel et al.l ()2005f ) and lSchoier et al] 
([2005). N2H+ has seven hyperfine lines that make up its 
(1-0) transition. We model the isolated 101-012 hyperfine 
line, which has a rest frequency of 93.176 GHz. The crit- 
ical density for LTE is estimated to be ricrit = 

1.4 X 10^ 

cm~^, using the relation Ucrit = ^ui/Kui where Aui is 
the Einstein A coefficient and Kui is the collisional rate 
(at an assumed kinetic temperature of 20K). As the iso- 
lated hyperfine line of N2H+ is optically thin it contains 
information about all the velocities along the line of sight 
and no asymmetric infall signature is expected. Nonethe- 
less, this line is a useful reference when compared to op- 
tically thick lines to separate self-absorption effects from 
intrinsic variation in the velocity profiles. 

We c alculate the a bundance of N2H+ using the re- 
sults of lAikawa et al.l (|200 5) who find an abundance of 
^N2H+ = 10~^^ relative to the H2 number density of 
the gas. The abundance of N2H+ has been proposed to 
rise when CO is froz en onto dust gra i ns as CO destroys 
N2H+. For instance |j0rgensen et al.l ()2004f ) found that 
as CO becomes depleted onto dust grains the abundance 
of N2H+ increases to A = 10~^. However, when we used 
this value for the N2H+ abundance above the CO freeze 
out density of n = 3 x 10"^ cm"^ (Bergin & Tafalla 2007), 
we found that the N2H+ emission became too bright and 
optically thick com pared to its observed behaviour. In 
lAikawa et al.l ()2005l ) the N2H+ abundances in collapsing 
Bonnor-Ebert spheres showed variations of at most a fac- 
tor of a few, and a slight decrease in abundance at the 



core centre. Therefore we will use the simple constant 
abundance model and leave a more detailed analysis of 
the N2H+ chemistry to future works. 

The second molecular tracer is CS, which is optically 
thick and traces the more diffuse gas in the envelope 
around the core. Molecular collisional ra tes and Ein- 
stein c oefficients for CS were calculated by iTurner et al.l 
()1992f ). We model the 2-1 transition of CS, which has 
a rest frequency of 98.981 GHz. The critical density of 
CS is ricrit = 3.2 x 10^ cm~^ at a collisional tempera- 
ture of 20K. As an optically thick line, CS exhibits self- 
absorption and is capable of producing a blue asymmet- 
ric collapse profile. The classic work on infall sign atures 
by [Zho3 ([1992) used CS. However, ITafalla et all (|2002h 
showed that CS is depleted at the centre of cores due to 
its freezing onto dust grains. They found the depletion 
could be described by an abundance prescription 

Acs(r)=Aoe-^^^)/^- (1) 

where Acs is the abundance, and n is the H2 number 
density. We adopt values of n^^ = 2 x lO^cm"^ for the 
depletion density and Aq = 4 x 10~^ for the low density 
abundance limit, b ased on the typ ical values given in 
ITafalla et ahl (j2002t ). Stahler fc YenI (|20i0 ) recently car- 
ried out line modelling of diffuse cores with and without 
CS depletion, and found that the models with depletion 
were a better match to observations, and so we adopt 
it here. Due to depletion the number density of CS is 
almost negligible in the densest regions of our cores. 

The third molecular species we consider is HCN, which 
is optically thick and traces the dense gas. The col- 
lisional rates for this spe cies have been calculated b y 
iGreen fc T haddeus (1974) and ' Dumouchel et al.l (j2010D . 
HCN is particularly useful because it has three well- 
separated hyperfine emission lines within its 1-0 level and 
each one has slightly different optical depths, with F(O-l) 
being the optically thinnest, F(2-l) the opti cally thick- 
est, an d F(l-l) intermediate in optical depth. iSohn et aJ] 
(j2007[ ) have recently pointed out the utility of this species 
as a diagnostic tracer of infall motions in dense cores 
and as such it is the key tracer we consider here for 
our filamentary cores. We use a constant abundance of 
^HCN = 3 X 10~^ for the molecular abundance of HCN, 
following iLee et ahl ()2004f ). Note t hat some v ariation in 
abundance is seen in core centres bv lLee et al.l (j2004) but 
for simplicity we neglect this here. The rest frequency of 
the central HCN hyperfine line is 88.632 GHz. 

To check that the cores are visible through the fila- 
ment, we first calculate the dust emission in RADMC- 
3D, using a dust-to-gas ratio of 0.01 and a constant dust 
temperature of 20K. Figure [2] shows that the embedded 
cores can be clearly seen in 850 fim dust emission. 

The line models are calculated for each molecular 
species at 256 evenly-spaced wavelengths, centred on the 
rest wavelength of the relevant emission line. The bot- 
tom panels of Figure [2] show the emission from Filament 
A at the rest wavelength for each of the three species 
(for HCN we show the F(2-l) hyperfine line). The N2H+ 
emission is mainly from the dense gas traced by the dust 
emission, whereas the CS shows more extended emission 
from the surrounding filament. The HCN emission is 
bright in both the cores and the surrounding regions. 
The second core, to the left of the central core, actually 
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TABLE 1 

The properties of the molecular tracers. The critical density for LTE is estimated using the relation = ^uil^ui 

WHERE Ay,i IS THE EiNSTEIN A COEFFICIENT AND Ky,i IS THE COLLISIONAL RATE COEFFICIENT AT AN ASSUMED KINETIC TEMPERATURE OF 

20K. 



Tacer Line Critical density Optically Abundance Law 

[n/njjj 

N2H+ 1-0 1.4 xlO^ thin 10-^° 

CS 2-1 3.2 xlO^ thick 4 x \{)-'^ e-'^nil (n^ = 4 X 10^ cm"^) 

HCN 1-0 2.6x10^ thick 3 x 10"^ 




line of sight [pc] 

Fig. 3. — The number density of CS (dotted line) and the gas 
number density (solid line) measured along a line through the core 
centre. There is a sharp rise in density at the location of the 
collapsing core (at a number density of around 3 x 10^ cm"*^). 
However at these densities the CS is frozen out and so it only 
traces the surrounding filament. 

appears brighter since it is more evolved and has a higher 
temperature. 

The abundance of CS is lower in the core centre than 
in the surrounding filament due to the effects of freeze 
out (outlined in the previous Section) . Figure [3] shows 
the number density of CS and the gas density along a line 
through the core centre. The densest regions of the core 
are effectively not contributing to the CS emission. The 
HCN and N2H+, however, follow the gas distribution. 
The gas densities are extremely low outside the filament, 
but increases steadily in the filament, and rises sharply 
at the center of the embedded core. As the density of the 
gas is so low at the edges of our box, we can be confident 
that we are integrating over the full length scale from 
which there are contributions to the total emission. 

4. RESULTS 
4.1. Line Profiles 

In Figure m we show the line profiles of the central core 
of Filament A, calculated for a beam with a FWHM of 
O.Olpc. The beam passes directly through the centre of 
the core and is viewed at a model inclination and rotation 
of zero. The dotted lines show the line profile of N2H+. 
As N2H+ is optically thin, all the emission is visible and 
the line profiles are generally Gaussian. This confirms 
that deviations from this profile in the other lines are 
due to optical depth effects. 

In Figure m the CS profile shows a slight red asymmetry 
rather than a blue one as the blue emission has been 
largely self-absorbed. However, higher intensity emission 
is seen from the blue side of the HCN F(2-l) line profile. 
The thinner F(O-l) line shows only a shoulder on the red 
side of a central peak, whereas the F(2-l) hyperfine line 
gives a clearer blue asymmetric profile. 



Generally the lines have more emission in their wings 
when observed in CS and HCN than in N2H+ emission. 
The bottom panel of Figure [2] shows the emission in each 
species from this filament. Due to their higher abun- 
dances, the emission from CS and HCN extends beyond 
the densest parts of the filament. This greater contri- 
bution of material from the outer edges of the filament 
increases the apparent line width of the core. The fila- 
ment velocities will be discussed in more detail in Section 

m 

Figure |4] shows a single viewing angle corresponding 
to an inclination and rotation of zero in our co-ordinate 
system. However, in such a filamentary and complex 
environment, the gas properties are likely to vary along 
each sightline. This is examined in Figures [5] and [S] which 
shows the HCN F(2-l) hyperfine line at various viewing 
angles around the filaments; the CS(2-1) lines are shown 
in Figures [7] and [8l In the left hand panels of the fig- 
ures the model is kept a the default rotation, ^=0°, and 
we view the model at 45° intervals in inclination. In 
the right hand panels, we keep a constant inclination of 
inc=90° and view the model at 45° intervals in rotation. 
We sample a total of 14 unique Hues of sight through 
each core centre. 

It is immediately apparent that there is a large de- 
gree of variability among the core line profiles depending 
on the viewing angle. Importantly, this is also true of 
sight-lines which pass through the same material but in 
the opposite direction. For example, the HCN lines at 
inc=0° and inc=180° in Filament A pass through the 
same gaseous regions but in the opposite direction, yet 
the former has a blue asymmetric profile and the latter 
a red. Even more striking examples in Filament A are 
the phi=90° and inc=270° sight-lines which do not even 
share the same general shape. This variability in the line 
profiles has immediate unfortunate implications for the 
interpretation of the optically thick line profiles for cores 
embedded within filaments. If a simple dynamical model 
(for instance that of a collapsing Bonnor-Ebert sphere) 
were used to predict observational line profiles, very dif- 
ferent interpretations of the dynamical state of this core 
would be deduced, purely due to line-of-sight effects. It 
is worth stating again that the core in the centre of each 
filament is the same regardless of viewing angle. 

Figure [5] contains sight-Hnes which show the expected 
blue asymmetry, but the majority of cases have no clear 
asymmetry at all. There are even a few cases that display 
a red asymmetry. A red profile is commonly taken as a 
signature of core expansion, or oscillation (e.g lKeto et al.l 
2006) as the reversed velocity field means that the red 
portion of the line profile originates from the higher den- 
sity part of a spherical core. Table [2] shows the total 
number of visible blue and red asymmetries in the HCN 
F(2-l) for our three embedded cores. We classify the 
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Fig. 4. — The emission line profiles from the embedded core at the centre of Filament A calculated using a simulated 0.0 Ipc half- width 
beam through the centre of the core. The dotted lines show the N2H+(l-0) optically thin line emission from the core. Higher intensity 
emission is seen from the blue side of the HCN F(2-l) line profile when the model is seen from this viewing angle. 



TABLE 2 

A SUMMARY OF THE CLASSIFICATION TYPES ASSIGNED TO THE 
HCN F(2-l) LINES FROM FILAMENTS A, B AND C USING THE LINE 
PROFILE SHAPES. DESPITE THE FACT THAT THE EMBEDDED CORES 
ARE COLLAPSING, A BLUE ASYMMETRIC LINE PROFILE IS SEEN IN 
ONLY 36% OF CASES. 



Filament 


Blue 


Red 


Ambiguous 


A 


5/14 


4/14 


5/14 


B 


7/14 


2/14 


5/14 


C 


3/14 


1/14 


10/14 


Total 


15/42 


7/42 


20/42 




36% 


17% 


47% 



cores as blue when there was a clear peak in the blue 
side of the profile and a clear lower intensity peak or 
shoulder in the red. When the more intense peak was on 
the red side of the profile, we classify the core as red. If 
both peaks were of a similar intensity, two clear peaks 
could not be distinguished within the line, or there was 
only one clear peak, the line profile is classified as being 
ambiguous. 

For HCN under half of the line profiles are blue, with 
a total of 35% of the sight-lines showing this asymmetry. 
Of the remainder, 17% of sight-lines show a red asym- 
metry, and 48% of cases show no clear asymmetry. For 
CS the interpretation is even more complicated as the re- 
sulting line profile is mainly determined by the dynamics 
of the filamentary envelope surrounding the core (since 
molecular freeze out of CS onto dust grains leads to very 
low abundances at the core centre, as shown in Figure 
[3]). Figure [71 shows that the CS line profiles are as likely 
to show a red asymmetric profile as a blue asymmetry 
and that the majority of the lines are hard to in terpret, 
particularly in filaments A and C. iPark fc Hong (fl998 ) 
showed that underlying 'dumpiness' within a core can 
flatten the resulting line profiles and superpose small ve- 
locity features onto the line where the maximum clump 
velocities lie. This complicates the interpretation of the 
line profile in irregular density fields, such as those in our 
models, and in some extreme cases can even lead to false 
infall expansion signatures. In our models we find only 
six line profiles which clearly resemble classical infall or 
expansion profiles in the CS sample, three (7%) red, and 
three (7%) blue. 

An alternative method of classifying cores is to use the 
normalised velocity difference SV between the optically 
thick and thin lines (Mardones et al..l997) . This is cal- 



culated by the equation 

SV = {VtMck - ytUn)imUn (2) 

where Al^^in is the full width half maximum (FWHM) 
of the optically thin line. We calculate this value for the 
isolated HCN hyperfine line and for the CS(2-1) lines by 
fitting a Gaussian profile to the N2H+ line to obtain the 
line centre and FWHM. We then find on which side of 
the N2H+ line centre the optically thick tracer peaks and 
fit a Gaussian to this component to obtain the optically 
thick component. Figure [9] shows the results of these 
calculations. 

A negative value of W is blue shifted and a positive 
value red shifted. An observation is considered to be 
a detection of true infall or expansion if \SV\ is greater 
than 5cr(JV). For HCN, 20/42 (48%) of the line profiles 
are classified as blue shifted collapse profiles, and 13/42 
(31%) of the line profiles are red expansion profiles. This 
is slightly higher than in the above classification by pro- 
file shape as some of the profiles previously classified as 
ambiguous are now classified as collapsing or expanding. 
It is interesting to note in Figure [9] that there is a gap 
in the distribution around the zero point. This could be 
either due to statistical effects or that our sample is bi- 
ased as it only contains dynamical collapsing cores and 
not any static cores. 

The W distribution for CS(2-1) is very slightly skewed 
towards the blue side. For CS there are 16/42 (38%) pro- 
files where W is indicative of collapse and 14/42 (33%) 
where bV is indicative of expansion. In addition to the 
high false positive rate of red profiles, many of the Gaus- 
sian fits to the CS(2-1) lines are too poor for the core 
to be classified. As in the case of the above classifica- 
tion by profile shape, the CS(2-1) line is a poor indicator 
of collapse for cores embedded in dense dynamical fila- 
ments. Given these detection statistics, it is clear that 
the velocities that the CS traces are largely unrelated to 
those within the dense core. Subsequently, we shall fo- 
cus on the HCN lines. As such we are taking the optimal 
case for detecting a blue asymmetric line profile from a 
collapsing core. 

All of the filaments contain a collapsing core, yet in 
more than 50% of cases, in even the most sensitive tracer 
(HCN), this collapse does not reveal itself as a blue asym- 
metry in the line profile. This leads us to the first con- 
clusion of this study: filaments can obscure signatures of 
collapse from their embedded cores. 
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Fig. 5. — The HCN F(2-l) line profiles of Filaments A (top), B (bottom) at different viewing angles. The central colour image of each 
panel shows the dust density in the plane in which the sight-lines pass through the core. The background grey-scale images show the large 
scale (0.4 pc diameter) dust emission map of the filament when viewed at the labeled angle. The position at which the outer images touch 
the central image denotes the orientation of the sightline. The line profiles of the central core at each angle are over-plotted on top of the 
grey scale image. A coloured dot shows how the profile was classified. The line profiles are calculated for a O.Olpc beam centred directly 
on the embedded core. In the left hand panels the rotation angle has a constant value of ^ = 0°, and in the right panels the inclination 
has a constant value of inc=90° 



4.2. Velocity structure of the filament 

Given the importance of the dynamical state of the 
envelope on the resulting core line profile, we examine 
in more detail the dynamics of the filaments. Figure 
[TOl shows the velocity field in three slices that intersect 
the core in Filament A, and a zoom into the centre of 
the x-y plane. Filament A is a turbulent sheet with a 
disordered velocity field within the x-y plane. There is 
no systematic velocity along the filament's major axis, 
and there are several regions that are contracting. In the 
y-z and x-z planes the diffuse material surrounding the 
filament has a clearer systematic velocity gradient and a 
larger absolute velocity. 

Sheets and filaments are well known to be one of 



the main sites where turbulent energy is transferred 
from large scales to small. For exam ple. [ Boldvrev et al.l 
(|2QQ2r ) extended the iShe fc Levequ3 (pJ9^ model of in- 
compressible turbulence to the supersonic compressible 
regime and found that making the assumption that dis- 
sipation occurred in micro-turbulent sheets achieved a 
good a greement with numer ical simulations. More re- 
cently, iSchmidt et al.l (|2QQ8l ) analysed high resolution 
turbulent boxes with various forcing scales and found 
that the most dissipative structures are intermediate be- 
tween filamentary and sheet-like structures. Therefore 
it is to be expected that shocks will dissipate energy 
and transform turbulence from high amplitude large scale 
motions in the diffuse gas into many low amplitude small 
scale motions within the denser filament. The velocities 
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Fig. 6. The HCN F(2-l) line profiles of Filament C at different 

shown in Figure [10] are therefore a natural consequence 
of filament formation from a turbulent molecular cloud. 

The zoom in to the central core of Filament A shown 
in Figure [TOb shows the detailed velocity structure sur- 
rounding the central embedded core. In this case the flow 
onto the central core is not radially symmetric, but in- 
stead there are twisting streams of material flowing onto 
the core. Behind the core there is even a stagnation point 
in the veloci ty field (sini ilar effects have been seen in 
iBalles teros- Paredes et al.] [2QQ3 and Klessen et al. 2005). 
Any sight line which passes through this stagnation point 
will diflFer substantially from the theoretical velocity pro- 
file shown in Figure [TJ 

Figure [11] shows a slice through Filaments B and C 
in a similar manner to Figure [10] In Filament B a fast 
moving flow of gas is pushing against a denser region of 
gas with lower mean velocities. The dense filament is 
found at the interface of the two regions and along this 
boundary turbulent whorls of gas can be seen. Along this 
interface material from the fast-moving diffuse flow will 
shock onto the surface of the filament. Filament C is once 
again formed at the interface between two flows. In this 
case, however, the contrast is even more extreme as the 
flows are converging. Both flows are supersonic which 
will lead to a rapid accumulation of dense material at 
their interface. 

Clearly the velocities in the filaments within which the 
cores are embedded are very different from the quies- 
cent enveloped envisage in Section [2] Moreover, these 
velocity flows are an intrinsic part of the filament as 
they result as a byproduct of its formation. Indeed it 
has been proposed that gravitational contraction and ac- 
cretion ar e drivers of turbulence on virtually all scales 
(iKl essen & Hennebell^ 120101). Many previous discus- 
sions of filamentary geometries in the literature have fo- 
cussed on either hydrostatic filaments ('Ostriker"1964') or 
their subsequent self- similar collapse (Hennebelle 2003; 
iTillev fc Pudri tz"^2003). The filaments in our large CMC 
simulations are dynamic objects that form as a result 
of turbulence and the associated shocks, and therefore 




angles. As in Figure [5] 



have a much mor e disordered velocity field (see also 
"Klessen et alJ|2000D. 

The effect that these disordered velocities in the plane 
of the filament have on the line profiles is illustrated in 
Figure [121 indicating the velocities and densities of the 
material along the inc=0°, phi=0° line of sight (this cor- 
responds to the -z axis shown in Figures [TO] and [TT]) . This 
figure provides examples of how the physical conditions 
in Figures [10] and [11] gave rise to the line profiles shown 
in Figure [5] 

The left panel of Figure [12] shows Filament A. Along 
this line of sight the stagnation point is behind the core, 
so there are no blue-shifted velocities. In other words, 
most of the mass the core accretion occurs from one di- 
rection. Consequently, when viewed from the forwards 
direction (inc=0°), the densest visible peak in the veloc- 
ity profile is on the blue side and when viewed from the 
the reverse direction (inc=180°) the densest visible peak 
is on the red side. This gives rise to a profile which is 
either red or blue asymmetric depending on the viewing 
angle. 

The central panel of Figure [12] shows Filament B. Here 
the velocity profile in the core is similar to that observed 
in Figure [TJ However, the region as a whole exhibits 
a large scale velocity gradient as the filament is swept 
along with a larger external gas flow. On the forward 
(inc=0°) side, the optically thick emission comes from 
the densest point of the blue-shifted velocities, and a blue 
asymmetric profile is produced. However on the reverse 
side there is optically thick emission at the same velocity 
in front of the core. This obscures the blue-shifted core 
emission, and so no blue asymmetry is seen. 

The right panel of Figure [12] shows Filament C. In this 
filament the velocity profile of the collapsing core joins 
smoothly on to the large scale convergent motions that 
are forming the filament. Here the converging gas flows 
first bring the filament together, and then at higher den- 
sities gravity leads to the formation of a dense core (see 
also Hennebelle & Chabrier 2008; Banerjee et al. 20d3)- 
When viewed from the reverse direction (inc=180°) the 
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gas distribution in front of the core does not overlap in 
velocity space with the central core. Consequently all 
the optically thick emission from the core is visible to 
the observer, and there is neither a self absorption dip at 
the centre of the line profile, nor any asymmetry between 
the blue and red side of the line. 

4.3. Line-widths 

Another interesting diagnostic is the linewidths of the 
N2H+ profiles. Table [3] shows the mean velocity disper- 
sion, (j{v)^ of the N2H+ lines when fitted by the following 
Gaussian function 



G = Ao exp 



(3) 



where Aq is the peak line brightness and vq the rest ve- 
locity. Table [3] also lists the velocity dispersion along the 
beam directly from the model. We calculate both the 
volume weighted and density weighted velocity disper- 



sions. The volume weighted dispersion is larger than the 
density weighted value, as expected from the fact that 
the core forms at the stagnation point of a convergent 
flow. To illustrate this. Figure [13] shows a histogram of 
the simulation velocities at equally spaced intervals along 
the beam for the filaments. In the case of Filaments B 
and C the velocity distribution is non- Gaussian and ap- 
pears bimodal, which reflects the converging flows of gas 
around the filament. In Figure [13] we also plot two sub- 
sets of the data; points with densities above 10^ cm~^, 
and points above 10^ cm~^. Both subsets occupy a nar- 
rower velocity range than the parent distribution, and 
the velocity range of the 10^ cm~^gas which is associ- 
ated with the cores typically is around the sonic scale. 
This explains why the 'observed' N2II+ linewidths of the 
cores are roughly sonic. 

That the velocity dispersions found from the N2H+ 
lines are typically coherent is in good agreement with 
recent observational studies of star- forming filaments. 
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Fig. 9. — Top The normalised velocity difference between the 
HCN F(2-l) hne and isolated N2H+(l-0) hne. Bottom As above 
but for the CS(2-1) line. Both distributions, but particularly the 
HCN, are skewed to the blue side of the distribution indicative of 
infall motions. 



iPineda et al.l (|2Q11[ ) carried out NH3 observations of a 
filament in the Barnard 5 region of Perseus and found 
it to be internally coherent, with a sharp transition to 
supersonic motions at its boundary. IHacar fc Tafallal 
|2Q11) observed filaments with embedded cores in L1517. 
They found that both the embedded cores and the dense 
gas in the filaments had subsonic velocities when ob- 
served in C^^O, SO and N2H+. The filaments in our 
models are formed from supersonic shocks. Before the 
shock (on the outside of the filament) the velocities are 
supersonic and the gas is diffuse, after the shock, the 
velocities are subsonic and the gas is dense in good 
agreement with these observatio ns. In general dense 
cores are coherent objects (e.g. ' Goodman et al.l 119981 : 
iKirk et al...2QQ7.) . An important feature of this simula- 



TABLE 3 

The mean line- width of the N2H+ line compared to the 

VELOCITY DISPERSIONS IN THE SIMULATION. THE VELOCITY 
DISPERSION OF THE LINE a(v) IS FOUND FROM A GAUSSIAN FIT TO 

THE LINE PROFILE FROM EACH VIEWING ANGLE, WHICH IS THEN 
AVERAGED TO FIND THE MEAN. THE VELOCITY DISPERSION OF THE 
SIMULATION IS FOUND BY CALCULATING THE BEAM WEIGHTED 
MEAN VELOCITY AND ITS STANDARD DEVIATION ALONG EACH 
SIGHTLINE, AND THEN CALCULATING THE AVERAGE OVER ALL THE 
SIGHTLINES. AS WELL AS THE VOLUME WEIGHTED AVERAGE WE 
ALSO CALCULATE THE DENSITY WEIGHTED AVERAGE. THE 
LINEWIDTHS ARE ROUGHLY SONIC DESPITE THE VOLUME WEIGHTED 
VELOCITY DISPERSIONS BEING SUPERSONIC. 



Filament 


Mean a{v) 
Line Profile 


Volume Weighted 


Density weighted 


A 


0.28 ±0.07 


0.45 ±0.31 


0.10 ±0.08 


B 


0.20 ±0.02 


0.94 ±0.67 


0.39 ±0.36 


C 


0.20 ±0.05 


0.91 ±0.47 


0.50 ±0.24 



tion is that we produce coherent cores purely by hydro- 
dyn amics without having to co nsider magnetic fields (see 
also Ballesteros-Paredes et al.'2003; Klessen et al.ll2QQ5l ). 

The hnewidths of the optically thick species appear 
broader than the N2H+ linewidths. Both N2H+ and 
HCN have similar Einstein coefficients and high critical 
densities. Consequently both species are tracing simi- 
lar gas; therefore this effect cannot be attributed to the 
underlying velocity field, but must instead be caused by 
optical depth effects. As the HCN emission is optically 
thick, emission from the dense core is obscured by the 
surrounding gas, which reduces the intensity of the cen- 
tral peak. The line wings, however, are less obscured and 
therefore make up a greater part of the total emission 
than in the optically thin case, explaining why the opti- 
cally thick lines appear flatter and broader. In the case 
of CS freeze-out this effect is compounded in the dense 
core, which further increases the relative contribution of 
emission from higher velocity gas in the filamentary en- 
velope. 

4.4. The surrounding line profiles 

For a core to be interpreted as collapsing, its line pro- 
file should differ from its surroundings. Figure [M] shows 
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Fig. 10. — The velocity vectors of material in three perpendicular slices through the centre of the embedded core in our original model 
of Filament A before radiative transfer. The background greyscale shows the densities in the slice (which has a depth of 0.002 pc.) 



the HCN F(2-l) hne profiles of the regions surrounding 
the central cores embedded within Filaments A and C. 
The cores are shown at an inclination and rotation of 
zero where each filament has a blue asymmetric line pro- 
file at its core. In each case the central core emission line 
is brighter than its surroundings. Moreover, the asym- 
metric line profile becomes less distinct as the distance 
from the central core and filament increases. The line 
profiles of the embedded collapsing cores can therefore 
be distinguished from their surroundings. 

Along the filament the line profiles have a greater blue 
asymmetry than in the gas outside the filament. This is 
particularly clear in Filament C which forms from col- 
liding gas flows. As shown by iMvers et al.l ()1996f ). a 



plane of two converging gas fronts also gives a blue asym- 
metric profile. In our filament the converging gas sheets 
first bring the filament together and then as the densities 
increase gravity causes the filament to collapse further. 
This reinforces the point made in the previous Section 
that the velocities within the filament as a whole also 
contribute strongly to the line profiles. 

We also examined the line profiles as a function of time 
to test how transient our results were. For the case of 
Filament A shown in Figure [lH a blue asymmetric profile 
was observed for at least a period of - 2 X 10"^ yr. After 
this period we expect outfiows to start affecting our core, 
an effect which is not included in our models, and so we 
do not extend our analysis beyond this time period. 
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Fig. 11. — The velocity vectors in slices through the central core in our model of Filament B top and Filament C bottom. The background 
greyscale shows the densities in the slice. The filaments are formed at the interface of converging flows. 



5. INTERPRETATION OF LINE PROFILES 

5.1. Blue profiles 

First we consider the HCN line profiles classified as 
blue in Section 14.11 Table |4] shows the mean of the peak 
intensity of the HCN F(2-l) and N2H+(l-0) lines from 
the three embedded cores. In each of the cores there 
is a greater peak intensity of emission from the blue 
asymmetric line profiles compared to the red and am- 
biguous profiles. (The standard deviations of the calcu- 
lated means have similar values at the extremes of their 
ranges, but there is still a clear trend towards higher in- 
tensities). However this is not true Jn^ the N2H+ line 
profiles. A survey by iLee fc: MversI (|2Qllf ) also found an 



increase in N2H+ emission in cores with a blue asymme- 
try in contrast to that found here. We account for this 
difference by noting that our sample only c ontains dense 
collap sing cores, but the sample studied bv lLee fc MversI 
()2Qllf ) may also contain cores in a less evolved state that 
are not yet collapsing. The authors propose that cores 
that are less bright in N2H+ are less dense and hence 
less likely to collapse. As all the cores studied here are 
already dense and collapsing they all have relatively high 
N2H+ emission. 

The trend of brighter peak temperatures in the opti- 
cally thick blue profiles can be understood by considering 
the origin of the emitting gas. In the blue cores there is 
typically no diffuse gas in front of the core. Therefore the 
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Fig. 12. — The line of sight properties averaged over our Gaussian beam along the inc=0° phi=0° direction. The line is solid in the core 
region (taken as a density above 3 x IC^ cm~^) and dotted where the gas is more diffuse. The velocity and density components along the 
line of sight are more complex than the simple symmetrical model shown in Figure ^ 
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Fig. 13. — Histograms of the distribution of velocities which make up the Gaussian beam along the inc=0° phi=0° direction for the 
filaments. The velocities in the clear histogram are sampled at equally spaced intervals along the beam. The patterned histogram contains 
only those points within the beam with a density above 10^ cm~^, and the solid histogram contains only those points with a density above 
10^ cm~^. The underlying velocity distributions are non-Gaussian and reflect the fllaments formation from bulk gas flows. However, the 
densest material typically occupies a narrow, typically sonic, region within this velocity space. 



TABLE 4 

The mean peak brightness temperature (in K) of the lines 
classified as blue, red and ambiguous when observed in 
HON F(2-l) IN Table [2] We also calculate the standard 

DEVIATION OF THE MEAN (nOT SHOWN FOR THE RED CORE IN 

Filament C as there is only one measurement.) The cores 
classified as blue have a greater hon brightness 
temperature since they primarily contain emission from 
the dense core rather than the more diffuse envelope. 



Filament 


Species 


Blue 


Red 


Ambiguous 


A 


HCN 


5.92 ± 1.84 


4.54 ±1.16 


4.07 ± 1.53 


A 


N2H+ 


1.14 ±0.32 


1.42 ±0.41 


1.82 ±0.64 


B 


HCN 


4.22 ± 1.07 


2.32 ±0.51 


2.48 ±1.07 


B 


N2H+ 


1.09 ±0.22 


1.14 ±0.05 


1.22 ±0.20 


C 


HCN 


5.51 ±0.89 


2.85 


3.44 ± 1.63 


C 


N2H+ 


1.20 ±0.36 


1.94 


1.25 ±0.40 



optically thick emission traces the high density core and 
the maximum brightness temperature is higher than in 
the red and ambiguous cases where the lines only trace 



the lower density gas. However, in the optically thin 
species the emission from the core is always visible and 
the foreground material represents an extra contribution 
to the total brightness. Therefore for an optically thin 
species the brightness can actually be higher in the red 
and asymmetric cases. 

We suggest that this may be a useful test in surveys of 
embedded cores. If the peak brightness is systematically 
higher across the survey in cores with a blue asymmetric 
line profile, but these identified cores do not show the 
same trend in their N2H+ profile, it is possible that some 
of the cores without a blue asymmetry are also collapsing. 

5.2. Red profiles 

Next, let us consider the existence of red asymmetric 
line profiles in our sample of collapsing cores. The red 
HCN line profile seen from Filament A at [inc=180°, = 
0°] in Figure [5] is an interesting case as it is the result of a 
one sided accretion flow onto the core. As the accretion 
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leading to one sided accretion onto the dense region. 

Another effect that can cause a red profile is random 
outward motions in the gas surrounding the core. As 
previously discussed, the peak brightness temperature of 
such profiles should be lower as the emission is coming 
from more diffuse gas. However, when the red asymme- 
try originates from a one sided accretion flow the bright- 
ness of the line is almost equal to that of the blue profiles 
because the emission originates from dense core gas. 

Red asymmetric profiles are seen in many observations 
of dense cores, for instance in the Gregersen et al. (1993) 
survey of Class cores, three out out of twenty-seven 
cores had a red asymmetry when observed in HCO+. In 
isolated cores surrounded by a hot confining medium, 
this line pro file is thought to indicate expansion or os- 
cillation (e.g lLada et al.|[2QQ3f ). However, the analysis of 
Filament A has allowed us to add a new potential expla- 
nation to this list. In cores embedded within filaments 
a red asymmetric optically thick line profile may be due 
to one-sided accretion onto the core. 
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Fig. 14. — The HCN F(2-l) lines of the region surrounding the 
central core of Filament A (left) and C (right). The background 
image shows the dust emission from the filaments. 



is only from one direction, the blue peak is higher when 
viewed from [inc=0°, = 0°] and the red peak is higher 
when viewed from [inc=180°, = 0°]. 

A one sided accretion flo w results wh e n an elongated 
core, such as those seen in ISmith et al.l (|2QTT1 ). IS more 
bound on one side than the other. Such a situation arises 
when a dense region of a filament becomes gravitation- 
ally unstable and, due to there being no pressure support 
until the formation of a protostar, the original geometry 
of the dense gas is maintained. Gravitational focussing 
encourages collapse at the ed ges of non-spherical objects 
(jHartmann fc Burkertl[2QQ7l ) and so the densest point of 
the core is frequently offset from the centre. The densest 
side of the core collapses quickly to the centre of poten- 
tial, but the other side of the core is accreted more slowly 



5.3. Ambiguous profiles 

Lines with no clear asymmetry are the hardest to in- 
terpret. In these cases turbulent line- widths within the 
filament can obscure optically thick core emission at the 
same velocity, and the profile appears identical to what 
would be expected if the central core were in equilibrium. 
As we showed in Section 221 turbulent motions are ubiq- 
uitous within filamentary structures as a result of their 
formation from shocks and shear flows. Consequently, 
the filaments in which the cores are embedded have a 
non-negligible velocity dispersion which can hide signa- 
tures of collapse from the optically thick lines of the cen- 
tral embedded core. Accordingly, a greater proportion of 
cores embedded within filaments may be collapsing than 
observations initially suggest. 

In the line profiles examined here, in more than 50% 
of cases the collapsing cores do not produce a blue asym- 
metric profile due to obscuration by their parent filament. 
If this were true for cores embedded within filaments in 
general, then a factor of 2.5 more cores may be collapsing 
than would be indicated by their line profiles alone. Such 
a substantial increase would suggest that dense cores are 
more unstable than previously thought. 

As discussed in Section 14. If the CS line profiles are 
both harder to interpret and are less reliable than the 
HCN observations. The vast majority of the CS emis- 
sion in our models originates from the filament enve- 
lope, partly as a consequence of freeze-out effects. Due 
to the turbulent motions of the filament, the line pro- 
files just trace the part of the filament along the line of 
sight. Any observational analysis of optically thick line 
profiles thought to originate from filamentary envelopes 
must take into account that information is only being 
obtained for one sight-line and that another sight-line is 
likely to give different results. 

5.4. Comparison to Observations 

Table [5] outlines the numbers of blue and red line 
profiles found in surveys of low mass cores. None of 
the surveys find the majority of the observed cores to 
have a blue asymmetry, but instead most find a frac- 
tion of above 30%. Only .Gregersen Evans. (j2000f ) 
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clearly excl u ded r ed asymmetries from their sample, and 
iSohn et all ^007t ) found 22% of their cores had a red 
asymmetry in at least two line transitions. 

When comparing to the observations several factors 
must be considered. The species used vary between the 
surveys and we have no HCO^ measurements in this 
work. Each core in the surveys is by necessity observed 
from only a single viewing angle, and it is not stated in 
the observations whether any of the cores are embedded 
within filaments. However, Herschel has revealed that 
many star-forming regions contain filaments which were 
previously not cleanly detected due to 'chopping' filter- 
ing out background structure within the observed cloud. 
It is therefore likely that some fraction of the cores stud- 
ied by these authors reside in more complex filamentary 
environments than was previously realised. 

We also note that caution is in order when interpret- 
ing the comparison to observations as various selection 
effects can be hard to account for in the observations 
compared to our simulations. Many of the observational 
studies considered here are biased towards cores previ- 
ously suspected to have infall motions. However our the- 
oretical study is also biased as we only contain cores that 
are unambiguously known to be collapsing as they go on 
to form stars. (The aim of our paper being to find if a 
core is collapsing, how often this can be deduced from 
the line profile when the core is embedded in a filament 
within a star forming cluster.) Given this fact, biased 
observational studies may actually the fairest compari- 
son. 

The observations of iLee et aH ()1999f ) primarily used 
CS which we find to not be a reliable tracer of collapse 
within our embedded cores. These authors found 29% of 
their cores to have blue line profiles, whereas our results 
were very ambiguous. There are several poss ible rea- 
sons for this discrepancy. Firstly not all of the Lee et al.l 
(|T999) cores are in filaments, some are more isolated 
cores, where the envelope of the core is i n deed collapsing 
and more easily detected. iLee fc MversI (|2Qlll ) revisited 
some of their sampled cores to examine their environ- 
ments and found that the majority (17/35) were associ- 
ated with filaments as considered here, a few cores were 
truly isolated (4/35) and the remainder (14/35) were 
surrounded by the cloud or at a cloud edge. Secondly 
our filaments are dense objects where CS is frozen out 
in the core, in less dense filament CS may be a better 
tracer. In observational studies it is frequently found 
that line profiles from different tracers vary widely (e.g. 
lShinnag"a in a similar manner to the theo- 

retical work pres ented here. 

The findings of' Sohn et al.l ()2QQ7f ) are the most directly 
comparable to this work as they focus on HCN observa- 
tions of dense cores. Their selected cores were nearby 
bright sources which were dense (n ~ 10^ — 10^ cm~^), 
compact (r ~ 0.05 — 0.35 pc) and had narrow N2H+ 
linewidths {6v ^ 0.2 — 0.4 kms~^). The surveyed cores, 
therefore, have similar phy sical properties to those stud- 
ied here. Sohn et aTl ()2007[ ) has the highest proportion of 
both red and blue line asymmetries of the surveys shown 
in Table [H In our models optical depth effects in the 
filament obscured the collapse signature from the cen- 
tral core in o ver 50% of cases. If the same effects were 
present in the iSohn et aTl ()2007f ) survey, it is within the 



bounds of possibility that the majority of their sampled 
cores were collapsing. 

Regardless of the species the root cause of the vari- 
ability with viewing angle is the intrinsic variation of the 
density and velocity fields with viewing angle. In our 
study the effect is most pronounced in CS emission as it 
was frozen out in the core, but there is also a high de- 
gree of variability within the HCN line profiles. Therefore 
we expect viewing angle effects to influence the observed 
line profiles of all tracer species, not just the ones stud- 
ied here. This leads to the possibility that all the surveys 
shown in Table [5] underestimated the number of collaps- 
ing cores, and dense cores may be more unstable objects 
than previously thought. 

6. CONCLUSIONS 

In this paper we consider three filaments with embed- 
ded collapsing cores from a large scale CMC simulation 
and perform radiative transfer calculations to model the 
line profiles. The filaments and their embedded cores 
form dynamically from bulk flows within the molecular 
cloud. The emission from each core is modelled using 
three species: optically thin N2H+(l-0), and optically 
thick CS(2-1) and HCN(l-O). Each model is observed 
from multiple viewing angles using a Gaussian beam 
through the core centre. Our main conclusions are as 
follows. 

1. The velocity fields within the filaments are disor- 
dered and contain turbulent velocities. The veloc- 
ities are imprinted on the filament as a result of 
its formation from turbulent converging flows and 
shocks. 

2. The optically thin N2H+(l-0) linewidths are 
roughly sonic, as shown in Table [3l 

3. The optically thick line profiles from the embed- 
ded cores are highly variable with viewing angle, 
as shown in Figures [5] and [71 

4. HCN lines, particularly HCN F(2-l) observations 
are more reliable indicators of core collapse than 
CS(2-1). In our simulated dense dynamical fila- 
ments CS(2-1) only traces the filamentary envelope 
around the core within which there is no systematic 
collapse in our models. 

5. More than 50% of cases show no blue- asymmetric 
collapse profile, due to emission from the turbulent 
filament surrounding the core. This effect could 
lead to an underestimation of the true number of 
collapsing cores in observational surveys. 

6. Sight-lines which trace the collapsing core are dis- 
tinguishable from sight-lines which mainly trace 
the filament, as they have a systematically higher 
HCN peak intensity than the rest of the sample. 
N2H+ emission does not show such a trend. 

7. In many sight-lines, a red asymmetric line profile 
arises from the embedded collapsing cores. This 
is either due to interference from the filamentary 
envelope, or one sided accretion onto the core. 

The data used for this work is publicly available by re- 
quest from the first author at the given correspondence 
e-mail address. 
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TABLE 5 

A COMPARISON OF SURVEYS OF INFALL IN LOW MASS CORES FROM THE LITERATURE. As A NOTE OF CAUTION NOT ALL SURVEYS CLASSIFY 
CORES AS 'blue' USING THE SAME METHOD. In THESE SURVEYS WE CLASSIFY CORES AS BLUE WHERE MULTIPLE COMPONENTS HAD A BLUE 

5V VALUE. When no comparable 5V values were given we used the number of infall candidates. (* Red profiles statistics 

WERE NOT DISCUSSED IN THIS PAPER.) 



Survey 


Species 


No. Cores 


Blue 


Red 


Gregersen et al. fl997) 


HCO+ 


23 


39% 


13% 


Greecersen & Evans r2000) 


HCO+ 


17 


35% 


0% 


Lee et aL C1999) 


cs 


69 


29% 


4% 


Mardones et aL (1997) 


H2CO,CS 


47 


32% 


* 


Andre et aL ^2007) 


CS, HCO+ 


25 


24-64% 


* 


_Sohn et al. (2007) 


HCN 


64 


43% 


22% 


This Work 


By shape 


HON 


42 


36% 


17% 


By 5V 


HCN 


42 


48% 


31% 


By 5V 


CS 


42 


38% 


33% 
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The line-of-sight densities and velocities are... 




In optically thick species only the nearest point at any 
velocity is seen. 



c 




self-absorption 




-V +v 



velocity 



Denser material emits at a greater intensity. 



